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ABSTRACT 

A model for pulsars with polar-cap magnetic flux density B antiparallel with spin f2 is 
described. It recognizes the significance of two elementary processes, proton production 
in electromagnetic showers and photoelectric transitions in ions accelerated through 
the blackbody radiation field, which must be present at the polar cap in the $7 • B < 
case, but not for pulsars of the opposite spin direction. The two populations are likely 
to be indistinguishable observationally until curvature radiation pair creation ceases 
to be possible. The model generates, and provides a physically realistic framework 
for, the polar-cap potential fluctuations and their time-scales that can produce mode- 
changes and nulls. The RRATs are then no more than an extreme case of the more 
commonly observed nulls. The model is also able to support the basic features of 
subpulse drift and to some extent the null-memory phenomenon that is associated 
with it. Unfortunately, it appears that the most important neutron-star parameter 
for quantitative predictive purposes is the whole-surface temperature T^, a quantity 
which is not readily observable at the neutron-star ages concerned. 
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1 INTRODUCTION 

The discovery within the last decade of further complex pul- 
sar phenomena, particularly the Rotating Radio Transients 
(RRATs; McLaughlin et al 2006) , has extended the problem 
of understanding these systems by introducing new time- 
scales. In a number of previous papers (Jones 2010a, 2011, 
2012a, 2012b; hereafter Papers I-IV) physical processes at 
the polar-cap surface in pulsars with spin ft and magnetic 
magnetic flux density B such that f2-B < have been exam- 
ined to see if they are relevant. The processes that introduce 
important time-scales are those associated with the forma- 
tion of electromagnetic showers by reverse-electrons from ei- 
ther electron-positron pair formation or from photoelectric 
transitions in accelerated ions. They are of no significance 
in n • B > pulsars with Goldreich- Julian charge density 
Paj < and outward electron acceleration. The work has 
assumed that the space-charge-limited flow (SCLF) bound- 
ary condition E • B = is satisfied on the polar-cap surface 
at all times. 

The aim of these papers has been to determine the com- 
position and energy distribution of the accelerated particle 
flux and to look for instabilities that might be relevant to the 
phenomena which are observed in the radio-frequency spec- 
trum. These properties can then be compared with those 



E-mail:p.jonesl@physics. ox. ac.uk 



needed to produce the observed emission spectra. This has 
a large bandwidth, of the order of 10^ Hz, and efficient con- 
version of kinetic energy to radio frequencies. For example, 
the peak power of the brightest pulses in PSR B0656-I-14 
(Weltevrede et al 2006a), is of the order of 7 x 10^'' erg s"\ 
which is equivalent to approximately 1Q^~^ MeV per unit 
charge accelerated at the polar caps. 

Growth of a collective mode able to transfer energy at 
this rate to the radiation field constrains the longitudinal 
effective mass, rrii-yf , of beam particles with mass rrii and 
Lorentz factor 7^. For a secondary electron-positron pair 
plasma, it is well known that this requires a low energy, 
7e ~ 100, which is also of the same order as that needed if co- 
herent curvature radiation, in a dipole-field, were the source 
of the observed radiation. Papers III and IV showed that, 
under the SCLF boundary condition, the creation of a re- 
verse fiux of electrons by photoelectric transitions in the ac- 
celerated ions limits the acceleration potential, analogously 
with the effect of electron-positron pair creation. Thus the 
outward particle fiux has two principal components: pro- 
tons formed in the electromagnetic showers and ions with 
Lorentz factors 7^ ~ 27a,z which are relativistic, but not 
ultra-relativistic as they would be in the absence of pho- 
toelectric transitions. They can have longitudinal effective 
masses that allow the rapid growth of the quasi-longitudinal 
Langmuir mode considered by Asseo, Pelletier & Sol (1990). 
As noted by Asseo et al, the quasi-longitudinal mode couples 
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directly with the radiation field and so introduces, in prin- 
ciple, a second source of coherent radio-frequency emission 
which is not present in fl • B > pulsars. 

Electromagnetic shower theory (see Landau & Rumer 
1938, Nordheim & Hebb 1939) makes it possible to calculate 
the total photon track length per unit interval of photon 
frequency. It is an almost linear function of primary elec- 
tron energy in the ultra-relativistic limit. Then known par- 
tial cross-sections for the formation and decay of the giant 
dipole resonance enable Wp, the number of protons formed 
per unit incident electron energy, to be estimated with ad- 
equate reliability for the present work. Proton formation is 
concentrated at shower depths of the order of 10 radiation 
lengths and it was shown in Paper I that the time Tp for dif- 
fusion to the top of the polar-cap atmosphere introduces the 
likelihood of instability in the composition of the accelerated 
particle flux. Estimates in Paper II, under more appropriate 
assumptions about the depth of the atmosphere, were that 
Tp would be of the order of 10^^ — lO" s. With this result 
in mind, it was ovidcut that the photoelectric transitions 
considered in Papers III and IV would introduce large-scale 
fluctuations in the the acceleration potential above the po- 
lar cap. The time-scale Tp is sufficiently long that the elec- 
tric field is simply given by solving Poisson's equation for 
the charge densities instantaneously present above the po- 
lar cap. It follows that the conditions necessary for coherent 
radio emission can appear or disappear with time-scales re- 
lated to Tp. 

Before proceeding further, it is appropriate to place this 
paper in the context of recent studies of the pulsar mag- 
notosphoro based on computational techniques in numeri- 
cal plasma kinetics. It is true that the inductance per unit 
length of the tube of open magnetic flux lines varies little 
with altitude above the polar cap so that the magnetosphere 
as a whole limits current-density time derivatives. The fur- 
ther hypothesis that the actual current densities arc deter- 
mined by the whole magnetosphere rather than by consid- 
erations specific to the small region immediately above the 
polar cap appeared in a paper by Mcstcl et al (1985) and 
was later extended by Beloborodov (2008). Solutions for the 
magnetic flux density external to a sphere, typically of radius 
^ 0.2Rlc, where Rlc is the light-cylinder radius, have been 
obtained by Kalapotharakos & Contopoulos (2009) and by 
Bai & Spitkovsky (2010) under the assumption of a force- 
free magnetosphere. Current densities can then be derived 
from the B-distribution as functions of basic parameters; the 
dipole moment, the rotation period P, and the angle ip sub- 
tended by fl and B. Those found by Bai & Spitkovsky devi- 
ate considerably from pojc, the Goldreich- Julian value, and 
were followed by one-dimensional time-dependent numerical 
studies of the polar cap with the E • B 7^ surface bound- 
ary condition (Timokhin 2010) and later with the SCLF 
boundary condition (Chen & Beloborodov 2013, Timokhin 
& Arons 2013) assuming fixed values of the time-averaged 
current density that are either larger than pojc or of the 
opposite sign. These produced evidence for the development 
of high potential differences and microsecond bursts of pair 
creation. 

All authors assumed the case fl ■ B > which is not 
the subject of the present paper. Some care will be needed in 
translating the results obtained by Chen & Beloborodov and 
by Timokhin & Arons to the f2-B < case. For this, the ba- 



sic unit of length in the description of the one-dimensional 
current and charge density, derived from the plasma fre- 
quency, is almost two orders of magnitude larger than for 
electrons. The one-dimensional potential must form a con- 
tinuous function with the three-dimensional potential de- 
rived from the Lense-Thirring effect (see Muslimov & Tsy- 
gan 1992). In this case, owing to the larger scale- length, 
there is no possibility of the particle backflow noted by Be- 
loborodov (2008) and Chen & Beloborodov (2013) so that 
only the stable zcro-tcmperaturc solution they refer to is re- 
alistic. The present paper does not accept the assumption 
of a precisely force-free magnetosphere on which the com- 
putational work is based. Instead, it retains the assumption 
that the current density is determined in the polar-cap re- 
gion and that it flows toward the light cylinder through a 
magnetosphere that is not precisely force-free. 

Paper IV attempted to describe, qualitatively, the rel- 
evance of physical processes at the surface of the polar cap 
to the complex phenomena that are observed in radio pul- 
sars, but a number of numerical estimates made earlier, in 
Papers I and II, have been superseded by later work. Thus 
the aim of the present paper is to introduce an elementary 
mathematical model of the polar cap which embodies those 
results obtained in Papers I- IV that remain relevant. The 
model itself and the information which has been derived 
from it are described in Section 3. Section 2 summarizes 
the physical properties of the polar-cap surface that are es- 
sential for its construction and also considers the non-static 
aspects of charge-to-mass ratio fractionation of atmospheric 
composition which were neglected in the earlier papers. 

It is the scale of the acceleration potential fluctuations 
which is the most significant feature displayed by the model. 
We relate this with the observed properties of the RRATs 
and nulls in Section 4 and with the phenomena of subpulse 
drift and null memory in Section 5. Our conclusions are sum- 
marized in Section 6 along with a brief discussion of observed 
phenomena on which, at the present time, our model appears 
to have no impact. 



2 ESSENTIAL FEATURES OF THE POLAR 
CAP 

It is assumed, unlike Papers I and II, that the actual polar- 
cap magnetic flux density is of the same order of magnitude 
as that inferred from the spin-down rate, which is usually 
less than the critical field m^c^/eh = 4.41 x lO" G. Thus the 
ion separation energy (Medin & Lai 2006) is small enough 
for the mass of the atmosphere at a polar-cap temperature 
Tpc ~ 10® K to be significant. 

This atmosphere is very compact and in local thermo- 
dynamic equilibrium (LTE); its scale height is of the order 
of 10~^ cm. The total mass is poorly known owing to its 
exponential dependence on the ion separation energy, but 
is possibly in an interval equivalent to 10^^ to 10^ radia- 
tion lengths. Thus it may contain the whole or some part 
of the electromagnetic showers formed by inward acceler- 
ated electrons. The extent of a shower is itself uncertain 
because the Landau-Pomeranchuk-Migdal effect is present 
at the energies, densities and magnetic fields concerned (see 
Jones 2010b). At depths immediately below the atmosphere, 
the state of matter is uncertain and could be either liquid 
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or solid (see Paper II). The nuclear charge is not known, 
but we assume the canonical value Z = 26. Electromagnetic 
showers reduce this to a mean Za and mass number A at 
the top of the atmosphere with an LTE ion charge Z. 



2.1 Atmospheric fractionation 

There is a fractionation of ion charge-to-mass ratio with the 
largest values at the top of the atmosphere. The relatively 
small number of protons created in showers are nowhere in 
equilibrium within the LTE ion atmosijhcrc and, under the 
influence of the small electric field E present, move out- 
ward and are either accelerated or, if their flux exceeds the 
Goldrcich- Julian current density pajc, form an atmosphere 
above the ions. Its scale height is 2kBT /rripg, at local tem- 
perature T and gravitational acceleration g. The chemical 
potential gradient that causes their motion within the ion 
sector of the atmosphere is initially mostly an entropy gra- 
dient, but at lower densities changes to that derived from 
the electric field. 

This fractionation is of central importance to the model 
described in this paper, but we have not previously consid- 
ered the adequacy of our elementary static treatment of the 
problem. Also, the presence of partial ionization means that 
the possibility of convective instability has to be considered 
(for a simple explanation see, for example, Rast 2001). Bear- 
ing in mind the functioning of a laboratory high- vacuum dif- 
fusion pump, it is also necessary to ask if the upward flux of 
protons is likely to carry with it sufficient ions to interfere 
with fractionation. 

In order to examine this, we require two transport re- 
laxation times. We employ approximate expressions using 
lowest-order zero-field perturbation theory, satisfactory here 
because the proton cyclotron energy quantum is small com- 
pared with the polar-cap thermal energy ksTpc. The first, 
for upward movement of a proton in the ion atmosphere at 
ion number density Nz is, 
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in which the function Fp is, 
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Here, Kop is the Debye wavenumber for the ions, 
«Lp = i^NzZ {Z + 2) 
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The equivalent ion relaxation time for movement relative to 
the proton atmosphere above the ions, if it exists, is also 
required. At proton number density Np it is, 
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in which Fz is given by equation (2) with the substitutions. 
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It is typically several orders of magnitude longer than r*"". 

The proton drift time from a depth z = with ion 
number density Nz to the top of the atmosphere is. 
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under the influence of a chemical potential gradient fixed by 
the ions. 
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equivalent to an upward directed force. The upward velocity 
is given by, 
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and the scale height of the atmosphere is, 
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(9) 



for A = 2Q,Z = Q,g = 2x 10^" cm s^^ and T = Tpc = 10'* K. 
Evaluation for Nz = 10^* cm'^^ gives I/t*'" = 2.5 x lO^'^ s"^ 
and a drift time of 1.0 s. This should be regarded only as a 
tentative order of magnitude because equations (1), (3) and 
(4) are strictly valid only at number densities much lower 
than 10^* cm~^. The position of the change of phase from 
gas to liquid or solid is also uncertain. 

The proton sector number densities are also sufficient 
for it to form an LTE atmosphere, but during phases of pro- 
ton emission into the magnetosphere it has the Goldreich- 
Julian outward flux. Thus as a flrst approximation, its ki- 
netic distribution function is isotropic in a frame of reference 
moving outward with a velocity pcjc/NpC until the proton 
number density Np becomes so small that the LTE condition 
breaks down. If this is not to interfere with fractionation, the 
force exerted on an ion by the upward flux of protons must 
be small compared with its chemical potential gradient in 
the static proton LTE atmosphere, which is. 



ArUpg — ZeE = ^ A — — \ rUpg. 
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A proton atmosphere able to maintain a typical Goldreich- 
Julian flux for one second would have a density at base of ~ 
10^^ cm~^ at which its ion transport relaxation time would 



be 1/rr 



5 X 10' 



"' giving a force several orders of 



magnitude smaller than equation (10) and of small effect on 
fractionation. 

The existence of convective instability at any density 
within the atmosphere is unlikely because the lateral motion 
implicit in a convective cell is strongly suppressed by mag- 
netic fields of the order of lO^'^ G. We refer to Miralles, Urpin 
& Van Riper (1997) for a full treatment of this problem. But 
even if convective cells exist, the presumption must be that 
the electric field E, reffecting the internal equilibrium of the 
adiabatically-moving volume, is still present within it. Thus 
the velocity of the protons, averaged over many circulations, 
would remain as calculated above and there would be no in- 
terference with fractionation. 



4 P. B. Jones 



2.2 Acceleration potential 



2.3 Photoelectric transitions 



The polar-cap radius, as in Papers III and IV, denotes the 
division between open and closed magnetic flux lines, and is 
that given by Harding & Muslimov (2001), 



Mo(0) = 



1/2 



(11) 



where P is the rotation period. Wc assume a neutron-star 
mass 1.4M© and radius R = 1.2 x 10® cm, for which /(I) = 
1.368. Our approximation for the electrostatic potential 3> 
is based on the Lense-Thirring effect described by Muslimov 
& Tsygan (1992), 



$(W, «) = TT {ul{z) - U^) {p{z) - PGj{z)) , 



(12) 



in cylindrical polar coordinates, at altitude z and radius 
u{z), for the specific case of a charge density p(z) that is 
independent of u and is only a slowly varying function of 
z. (It is almost identical with the potential that would be 
present given a time-independent outward flow of electrons 
under SCLF boundary conditions in fJ ■ B > pulsars.) 
It also assumes approximate forms, valid at altitudes well 
within the light-cylinder radius, for the more precise charge 
densities which were given by Harding & Muslimov (2001). 
These are independent of vi{z) as required by equation (12) 
and are. 
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in which V is the angle between f2 and B, k = 0.15 is the 
dimensionless Lensc-Tliirring factor, and rj = {1 -\- z/ R). 

Equation (12) is certainly a satisfactory approximation 
to the true potential at altitudes z^ uo, but for lower val- 
ues of z, we must assume that the H-M potential changes 
continuously to the one-dimensional potential that exists at 
« -C Wo, which was originally described by Michel (1974) 
and then further investigated by Mestel et al (1985) and 
by Beloborodov (2008). We have already noted, in Section 
1, that the length scale associated with the form of the 
one-dimensional potential is large, in the ■ B < case. 
Thus equation (12) is a fair representation of the true po- 
tential, which would remove any possibility of the backflow 
described by Beloborodov. 

Under the assumption that equation (12) is valid, the 
condition p(0) = pajiQ) ensures that the SCLF condition 
E • B = is satisfied at all u on the polar-cap surface. The 
maximum potential available for acceleration above the po- 
lar cap occurs if p[z) is independent of altitude as would be 
the case if the particles were ultra-relativistic and there were 
no charge-separating interactions. Expressed in convenient 
energy units, it is. 



Vmax{u,00) = 



2n^R^KeB 



1.25 X loM 1 - ^ 



^ GeV (15) 



The blackbody temperatures responsible for photoelectric 
transitions are those of the polar cap, Tpc, and of the whole 
surface, Tg. The polar cap is not significant at altitudes 
z » wo(0) because the photon Lorentz transformation to 
the ion rest frame becomes unfavourable. There is a little 
ionization at z < h ~ 0.05i? to a mean charge Zh but the 
ion Lorentz factors there are small, given the SCLF bound- 
ary conditions, as is the contribution to the total reverse- 
electron energy per ion. We assume a fixed value for this 
component of ~ 20 GeV. Lorentz transformations of 
whole-surface photons are much more favourable at higher 
altitudes and produce the major part of the total reverse- 
electron energy per ion accelerated, e = en + Ss- Ionization 
may or may not be complete and we define the mean final 
charge as Zoo- Wc have shown previously that the flux of 
shower-produced protons reaching the top of the neutron- 
star atmosphere at any point u is given by. 



f 

J — oc 



dt7p(t-t')i^(u,tV"(u,^'). (16) 



per unit charge on a flux line with radial coordinate u. 



in terms of the ion flux (see equation (20) of Paper IV). 
The first component cannot exceed the Goldreich- Julian 
flux; the remainder J*" accumulates at the top of the atmo- 
sphere. Paper I assumed the neutron-star atmosphere to be 
of negligible depth so that /p was assumed to be the stan- 
dard diffusion function. But for the depth considered in Par 
per IV and here, an expression based on the drift time given 
by equation (7) is a better approximation. If the proton at- 
mosphere is exhausted, ion emission commences extremely 
rapidly in order to satisfy the SCLF condition E • B = at 
the polar-cap surface. In general, the screening of any elec- 
tric fleld in the fl • B < case for which a positive charge 
density is needed occurs preferentially through ion or proton 
emission in a single relativistic-particlc transit time, whereas 
the process of electron-positron pair multiplication requires 
many transit times. 

Protons accelerated to 7p ~ 10'' have only a very small 
probability of creating electron-positron pairs through inter- 
action with blackbody photons and, with ions of the same 
energy per unit charge, would have only negligible growth 
rates for the quasi-longitudinal Langmuir mode. There are 
three ways in which particle beams capable of giving strong 
coherent radio-frequency emission might be produced. The 
two considered here in the first instance are as follows. 

(a) Vjnax is so large that self-sustaining curvature- 
radiation (OR) electron-positron pair production occurs. A 
permanent proton atmosphere exists over much of the po- 
lar cap giving a primary current density of protons and 
positrons which, at least for times long compared with uo /c, 
may be in a steady state. A plasma of low-energy secondary 
electrons and positrons forms. 

(b) The potential is so reduced from Vmax by photo- 
electron backflow that ion and proton Lorentz factors are 
either in a region that allows rapid growth of the quasi- 
longitudinal mode or are of magnitude such that downward 
fluctuations to the necessary values can occur. In this case, 
the Lorentz factors arc such that the mode wave- vector com- 
ponent perpendicular to B is unlikely to be negligible so that 
coherent radio-frequency emission need not be exactly par- 
allel with local flux lines. 
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The third set of processes are more obscure in the case 
of i7 • B < pulsars and concern inverse Compton scattering 
(ICS) of blackbody photons above the polar cap. Pair pro- 
duction by the conversion of outward-moving high-energy 
ICS photons is known to be significant in • B > pulsars 
(see Hibschman & Arons 2001, Harding & Muslimov 2002) 
even if dipole-field geometry is assumed. Also, deviations 
from such a field, to the extent that they exist, can greatly 
enhance pair densities (Harding & Muslimov 2011). But in 
• B < pulsars, high-energy ICS photons arc directed in- 
ward and the sources of outward-moving positrons to scatter 
blackbody photons are more obscure. Clearly, a fraction of 
the photons may have sufficient rnorncnturn transverse to B 
to produce pairs at low altitudes. A further possible source 
of low-altitude pairs is the conversion of neutron-capture 
7-rays at the top of the atmosphere. Showers produce ap- 
proximately the same numbers of neutrons and protons, but 
reliable calculation of the probability that a neutron would 
reach the top of the atmosphere is not easy. It is also not 
obvious how the density of electron-positron pair production 
via these processes would respond to fluctuations in the ac- 
celeration potential. However, although the model described 
in Section 3 has been considered in the context of cither low- 
energy ion-proton beams or CR pair production, it is neces- 
sary to bear the more obscure pair-production processes in 
mind. 



3 THE POLAR-CAP MODEL AND ITS 
PROPERTIES 

The previous Section summarized those specific physical 
properties of the polar cap that are important for the model, 
but we shall also give a brief description of its framework. 
This is addressed primarily to the typical radio pulsar with 
an age of the order of 1 Myr, or older, and not to the small-P 
case in which continuous pair production is anticipated. 

Photoelectric transitions in accelerated ions arc a source 
of electrons which flow inward to the polar-cap surface and 
partially screen the potential given by equation (12). Data 
presented in Table 1 of Paper IV, for specific values of B and 
Ts, gave Zoo, the final charge of an accelerated ion, and the 
energy of its reverse electrons €s , as functions of the acceler- 
ation potential V{u, oo) experienced by that ion. Then with 
the inclusion of the small energy eh arising from photoelec- 
tric transitions at altitudes z < h, as described in Section 
2.3, the total energy e = en + es gives the number of protons 
created in the electron showers. This is the approximately 
linear function eWp — KZoo (see Papers I and II). It is con- 
venient here to define K as the number of protons formed 
per unit ion charge. (In Papers I and II, it was defined as 
the number of protons per unit nuclear charge of the ion.) 
The model described here uses an elementary form of equa- 
tion (16) which describes the temporal distribution of the 
protons reaching the top of the LTE atmosphere. But the 
distribution of electron back-fiow and the values of Zoo over 
the polar-cap also define through solution of Poisson's equa- 
tion, at any instant, the actual electrostatic potential which, 
of course, is not identical with equation (12) because p is not 
independent of u. The model then finds by relaxation pro- 
cedures finite-element approximations to the self-consistent 
functions Zoo{u,t) and V{u,oo,t) for the polar cap. This is 



circular, with radius given by equation (11), which assump- 
tion is not essential but is made for ease of calculation. 



3.1 The model 

The polar cap is divided into ris = elements of equal area 
by first defining n annuli of outer radius iuo/n, where i = 
l....n and then further dividing each of those into 2i — 1 equal 
azimuthal elements. Within each annulus, the elements are 
displaced, by a random azimuthal angle, with respect to 
the axis </> = 0. It is assumed that the charge density is 
independent of u within any element. Equation (16) is much 
simplified by the model assumption fp{t) = 6{t — t' — Tp) 
giving. 



r{t) + rit) = K{t - Tp)r(t - t^), 



(17) 



within any element, so that the state of the element alter- 
nates between proton and ion emission. The total number 
of protons produced at time t on unit area of surface at the 
top of the atmosphere following an ion-emission interval of 
length Tp is. 



i: 



dt'K{t')r{t') 



(18) 



The ion current density is J'{t) = Nz{t)cZoo{t)e and the 
Goldreich- Julian charge density is pGj{h) = Nz{t){2Zh — 
Zoo{t))e. Therefore, if we define the length of the proton 

phase by Krp = Qp/pcjc, 



Ktp 



dt' 



K{t')Zoo{t') 

2Zh-Zoo{t'y 



(19) 



An event is defined as a change of state in any element from 

proton to ion emission or vicc-versa. The length of an ion 
phase is Tp, which is assumed constant. The length of the 
following proton phase is KTp, and is a function of the state 
of the whole polar cap within the integration time over the 
preceding ion phase. 

Table 1 of Paper IV gave values of and of Zoo as 
functions of a cut-off potential Vc for different values of the 
surface magnetic fiux density -B12, in units of 10^^ G, and the 
whole-surface temperature Ts . For reasons which we outline 
later, we assume a surface nuclear charge Zs = 10. The cut- 
off potential is a representation of the effect of photo-electric 
transitions on the maximum possible acceleration potential 
difference Vmax given by equation (15). From Table 1, with 
some additional values, we have been able to tabulate the 
functions (.s{V) and Zoo{V) for the interval < y(oo) < 
Vmax- They arc not strongly dependent on u and we can 
assume that they are functions only of the magnitude of 
y(u, 00), the potential at the upper end of the polax-cap 
acceleration zone which is represented here by a cut-off at 
?7 = 4. In the model, it is assumed to be the potential on the 
central axis of an clement. The excess charge density within 
an element during the ion phase is. 



2[Zoo — Zh) 
2Zh — Zoo 



PGj{h), 



(20) 



and is the source of a potential deviation downward from 
Knax(u, 00). In order to calculate this with some economy, 
we use the approximation on which equation (12) is based, 
specifically that at an altitude 2: >> uo, a section of the long 
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narrow open magnetosphere can be approximated locally by 
a cylinder of constant radius uo{z). Then the Green function 
satisfying the SCLF boundary condition for a line source at 
cylindrical polar coordinates (u, (f)) is, 

rCi. ..'1 -\r,( "° + ^^"'^ ~ ^^0""' cos((?!) - (j)') \ 
' " ^ - \yulv? + ulu''' - 2uluu' cos(^ - ■ 

The potential on the axis of an element derived from the 
charge excess within that element is the most important 
factor and is obtained by numerical integration using the 
Green function. The potential generated on its axis by any 
other ion-phase element is found directly from the Green 
function by using a line approximation for the excess charge. 
The potential within a proton-phase element is of no interest 
and is not calculated. 

After each event, the Zoo in all ion-phase elements, also 
referred to as ion-zones, are recalculated by a relaxation 
procedure so that there is consistency between the Poisson 
equation solution for V(u, oo) in terms of all the Z^o and 
the tabulated function Zoo{V) obtained from Table 1 of Pa- 
per IV which is based purely on the photoelectric transition 
rates. This also gives the function and hence K which, 
with the values of Zoo and their times of duration, then al- 
low calculation of the integral for K and the duration of the 
subsequent proton phase. 

The initial state of the system has all elements in the 
proton phase. All Us elements, selected in random order, 
are then assigned sequential times for transition to the ion 
phase. The interval between any two adjacent times is XTp, 
where a; is a random number in the interval < a; < Xmax- 
Our initial choice was Xmax = 0.5. Throughout the calcula- 
tion, a list is maintained, in temporal order, of the time of 
the next event in each of the na elements. This procedure 
has the advantage that the calculation of every ion phase 
is entirely self-contained so that currmlative errors are not 
carried forward. The model has been run for intervals equiv- 
alent to real times of the order of ten days. 



3.2 Model results 

For a given random initial state, the subsequent states of the 
model polar cap are not, of course, random but are chaotic, 
and variation of Xmax appears not to affect their character. 
The proton production parameter Wp is a slowly varying 
function of B but wo have chosen the conservative value 
Wp = 0.2 GeV~^ and an ion charge Zh = 6 as in Paper IV. 
The numbers of elements have been in the interval 10^ ^ 
ns ^ 4 X 10^ 

The usual procedure has been to let the model run for 
a time 10* before samjjling its state at intervals of Tp. The 
presence of the instability described in Paper I has been 
confirmed. The system has shown no sign of settling down 
to other than a chaotic state. Its main characteristics are as 
follows. 

(i) Fluctuations in the central potential V{0, oo) and 
in the number of ion-emission elements can be very large 
and are dependent, principally, on the whole-surface tem- 
perature Ts but also, to a lesser extent, on the parameter 
B\2P^^ which scales Vmax- 

(ii) The elements most likely to be in an ion-emission 
phase are those near the periphery uq. 



Table 1. The first two columns give values of the rotation pe- 
riod P and whole-surface temperature Ts . The next five columns 
histogram the a<;celeration potential V = V{0, oo)/Vmax{0, oo) 
sampled after 5000 successive intervals of Tp in order to show the 

scale of the fluctuations that arc present in the model for various 
P and Ts - The bin size is 0.2V and tiie number of elements here 
is rts = 100; the polar-cap magnetic flux density is B12 = 3.0. 
The final column gives Pvn, the probability that the whole polar 
cap is in the proton phase at any instant. 



p 








no) 






Pvn 


s 


10^ K 


0-1 


1-2 


2-3 


3-4 


4-5 




1.0 


1.0 


0.000 


0.004 


0.008 


0.021 


0.967 


0.000 




2.0 


0.011 


0.013 


0.019 


0.061 


0.784 


0.112 




4.0 


0.018 


0.064 


0.144 


0.302 


0.464 


0.008 


2.0 


1.0 


0.000 


0.000 


0.000 


0.000 


1.000 


0.000 




2.0 


0.000 


0.000 


0.001 


0.146 


0.854 


0.000 




4.0 


0.040 


0.247 


0.384 


0.280 


0.048 


0.001 


3.0 


1.0 


0.000 


0.000 


0.000 


0.000 


1.000 


0.000 




2.0 


0.000 


0.000 


0.004 


0.388 


0.608 


0.000 




4.0 


0.077 


0.464 


0.364 


0.085 


0.009 


0.000 



(iii) Peripheral ion-emission elements have a significant 
autocorrelation function in the angle 4>, unaffected by run- 
ning time, indicating the formation of clusters. 

At this stage, we should note that the model as precisely 
defined by equations (17) - (19) does not produce, at any 
point on the polar-cap surface, the mixture of protons and 
ions which is required for growth of the quasi-longitudinal 
Langmuir mode. But the growth rate is only a slowly varying 
function of the proton-ion ratio (see Papers III and IV) and 
we do not doubt that, given a physically realistic form of 
the diffusion function fp in equation (16), there will be a 
sufficient interval of time, of the order of Tp, within which 
the current density contains both components. 

The potential fluctuations given in the Table are, of 
course, dependent on Us. The autocorrelation function in 
indicates that a smaller of Us should be preferred and, within 
the limited framework of our model, would be a better ap- 
proximation to reality. But Ua = 100 is the smallest value 
for which our approximation for the interaction between two 
elements can be reasonably adequate. Table 1 therefore un- 
derestimates the true scale of fluctuation. The values of B12 
and P have been chosen to be representative of the typi- 
cal pulsar exhibiting nulls as listed by Wang, Manchester & 
Johnston (2007) . The scale of the downward fluctuations in 
V is, as might be anticipated, very strongly dependent on 
Ta and hence, presumably on pulsar age. At Ts = 10* K, 
significant downward fluctuations would be present only for 
very large values of Bi2P^^, the parameter that scales Vmax- 
With increasing rotation period, downward fluctuations also 
increase. 

The most significant feature of the model is the size 
of the potential fluctuations that occur for larger values of 
Ta. But this is superimposed on a further instability, with 
medium time-scale, described in Paper II. Showers reduce 
the atomic number from Z, possibly Z = 26, to Za at 
the top of the atmosphere. Provided we ignore capture of 
shower-produced neutrons and subsequent /?-decay of the 
neutron-rich nuclei formed, this liberates Z — Zs protons 
which diffuse with time-scale Tp. At any instant, there is 
naturally a distribution of Za values, but over a time long 
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compared with the ablation time (the interval in which one 
radiation length of matter is removed from the polar cap at 
the Goldreich-Julian current density) 

rH=2.1xl0^f Zl^^£±\ (22) 

\ZBi2 ln( 12^1/25- 

the average quantities must satisfy 

Z-{Zs) = {KZs). 

This means that Zs and K defined by equation (19) can- 
not be independent variables when considered over medium 
time-scales of the order of Tri- The reason is that large val- 
ues of K within some interval of time imply reduction of 
Z at the shower maximum to very small Zs so producing 
nuclei able to move to the top of the atmosphere either by 
diffusion on a time-scale possibly some orders of magnitude 
larger than Tp, or by Ray leigh- Taylor instability. In partic- 
ular, values smaller than ~ 5 have a high probability of 
being completely ionized in the LTE atmosphere and the 
consequent absence of a reverse-electron flux stops shower 
and proton formation. The condition of these surface layers 
is obviously not well understood, but in Paper II a case was 
advanced that it would be one of instability rather than a 
steady-state value of Zs- This is the case for our model as- 
sumption of a compromise value Zs = 10 and is also the 
basis for the observed phenomena, mode-changes and long 
nulls, with time-scales of the order of Tri- 



4 RRATS NULLS AND POTENTIAL 
FLUCTUATIONS 

Our proposal is that the potential fluctuations found in the 
model provide a basis for understanding why the extent 
of null lengths and fractions seen in both pulsars and the 
RRATs occur quite naturally. 

4.1 The RRATs 

We refer to Burke-Spolaor & Bailes (2010), Kcane et al 
(2011) and Keane & McLaughlin (2011) for reviews of the 
observational data on these sources. Values of fon, the frac- 
tion of periods in which a pulse is detected, are broadly in 
the interval 10~^ < fon < within which their distri- 

bution is approximately uniform in log /on. The bursts of 
emission are short, usually one but occasionally several pe- 
riods, and are consistent with being of the same order of 
magnitude as Tp. The RRATs listed in Table 3 of Keane 
et al are all quite distant with the exception of J1840-1419 
which also has a period P — 6.6 s and a value of the accel- 
eration potential parameter B\2P~^ = 0.15, that is, below 
the cut-off value 0.22 which can be found from the ATNF 
catalogue (Manchester et al 2005). 

Thus Weltevrede et al (2006a) proposed that the RRAT 
pulses are simply analogues of the giant pulses seen in the 
nearby pulsar B0656-I-14 and that their emission, apart from 
the giant pulses, is unobservable owing to distance. How- 
ever, there is a difference whose significance becomes obvi- 
ous in view of the model described here. With the exception 
of J1554-5209, the 14 RRAT Usted with surface magnetic 
fields by Keane et al are well below the threshold for CR 



secondary pair production as defined by Harding & Mus- 
limov (2002). EYom Fig.l of their paper, this condition can 
be represented approximately in terms of the parameter 
X = SiaP"^-^: the threshold is given by the critical value 
Xc = 6.5. PSR 0656+14 lies well above this threshold, which 
assumes dipole-field fiux-line curvature, as does the excep- 
tional RRAT J1554-5209, and it is reasonable to suppose 
that these pulsars are capable of supporting CR secondary 
pair creation over at least a fraction of their polar caps. 
Therefore, we suggest that the giant pulses of B0656+14 
are produced by process (a) through a rare upward fluctu- 
ation in V toward Vmax and its more general emission by 
process (b). Some RRATs may, of course, have sufficient 
non-dipolc flux- line curvature to enable process (a), but dis- 
counting this possibility, the RRATs listed by Keane et al 
are so far below the Harding-Muslimov threshold that we 
suggest their isolated pulses are a consequence of rare but 
large downward fluctuations in V. 

It might be argued that, owing to the boundary condi- 
tion satisfied by equation (12), Lorentz factors 7a, z small 
enough for growth of the quasi-longitudinal Langmuir mode 
always exist, in principle, at u close to wo so that there 
should be some coherent emission at all times. But in this 
case it is not obvious that the lateral (perpendicular to B 
and to u) depth of the particle beam would be suitably large 
in relation to the wavelength of the mode to allow the nec- 
essary growth rate. There is an obvious need for some ana- 
logue of the rclativistic Penrose condition (see Buschauer 
& Benford 1977) extended to allow for a beam with lateral 
velocity gradients and limited lateral depths. The question 
of the conditions necessary for observable emission will be 
further considered in Section 5.1 in relation to the structure 
and drift of subpulscs. 

The ages and surface magnetic fields of RRATs given 
by Keane et al can be compared with those of the nulling 
pulsars listed in Tables 1 and 2 of Wang, Manchester & 
Johnston (2007). It is immediately seen that RRAT fields 
are typically an order of magnitude larger than those given 
in the ATNF catalogue (Manchester et al 2005) for the Wang 
et al list. But the distributions of the acceleration param- 
eter Bi2P~^ are very similar because the RRAT have long 
periods. The age distributions are also similar though both 
are very wide. 

The relation between age and proper-frame surface tem- 
perature Ts is, unfortunately, obscure in the region assumed 
in Table 1, which is well within the photon-cooling epoch. 
Potekhin & Yakovlev (2001; Fig. 7) have shown that increas- 
ing surface fields accelerates cooling, but an even more im- 
portant factor would be the presence of law-Z elements in 
the outer crust, possibly from fall-back, which have the same 
effect (Potekhin et al 2003). 

4.2 Null lengths and surface temperature 

With reference to Table 1, it is possible to sec qualitatively 
how nulls change with age during the cooling from high-T^ 
to low-Ts. We consider pulsars which lie below the pair- 
creation threshold. Initially, values of V are usually small 
enough to give the mode growth rates necessary for emis- 
sion by process (b) (see Section 2.3). Upward fluctuations 
of V increase the ion and proton Lorentz factors and, ow- 
ing to the exponential dependence of amplitude growth on 
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them, can produce short nuUs. Upward fluctuations become 
more frequent as cooling proceeds, giving higher nuh frac- 
tions 1 — fon- Eventually, cooling reaches the stage at which 
the mean V approaches unity and downward fluctuations 
arc required to produce the growth rates necessary for ob- 
servable emission. Null fractions are then near unity and, 
as the limit is approached, there is no observable emission. 
We emphasize that these fluctuations are superimposed on 
medium time-scale fluctuations (see Section 3.2) in the sur- 
face atomic number Zs which arc unlikely to bo uniform over 
the whole polar-cap surface. Values smaller than ~ 5 pro- 
duce a negligible reverse-electron energy flux and hence, if 
present over a substantial area of the polar cap, much re- 
duced deviations of V from unity. It is proposed here that 
these are responsible for null and burst lengths more nearly 
of the order of th rather than Tp. 

Isolated Neutron Stars (INS) are a small group of radio- 
quiet thermal X-ray emitting sources positioned close to the 
RRAT in the P - P plane (see Keane ct al 2011) but at 
slightly greater periods. The six sources for which timing 
solutions exist have been listed by Zhu et al (2011). Given 
their kinematic ages, the observer-frame temperatures are 
large (~ 10^ K) compared with those considered in Table 1 
but, with the exception of the anomalous J0420-5022, their 
B\2P^^ are typically close to the general cut-off value of 0.22 
found from the ATNF catalogue (Manchester et al 2005). If 
they form part of the f2 • B > population, they should sup- 
port pair formation by the ICS mechanism; for fJ • B < 0, 
extrapolation from Table 1 shows that small V would be ex- 
pected satisfying the condition necessary for rapid growth 
of the ion-proton beam instability. However, the fact that 
no radio emission has been observed from this small num- 
ber of sources is, perhaps, not surprising in view of their 
anticipated small beaming fractions compared with the 47r- 
observability of the whole-surface X-rays. 



5 SUBPULSE DRIFT AND NULL MEMORY 

Subpulse drift, and the null memory which, in some pulsars, 
is observed with it, is obviously an important diagnostic of 
polar-cap physics. Unfortunately, our model is incomplete in 
that it does not have the capacity to spontaneously exhibit 
this phenomenon following the randomly-constructed initial 
state described in Section 3. Its flnite-element structure and 
the very elementary approximation represented by equation 
(17) are the reasons for this. However, the phenomenological 
model developed by Dcshpandc & Rankin (1999) from the 
classic Ruderman & Sutherland (1975) model has proved so 
useful that any physical polar-cap model must be able to 
support it. 

5.1 Nulls and null memory 

It was noted in Section 3 that the model system settles down 
to a chaotic state, but with ion-emission elements most likely 
to be near the periphery uo and an autocorrelation function 
in the azimuthal angle <i> indicating the formation of clus- 
ters. This is unsurprising, because ion-emission elements at 
smaller u have larger reverse-electron energy fluxes and so 
lead to longer phases of proton emission. Thus our model 
generates quite naturally the conal structure that is required 



for the Dcshpandc- Rankin carousel model and we can there- 
fore proceed to examine the extent to which it is capable of 
supporting its organized subpulse motion. 

For convenience, we shall consider a circular path of ra- 
dius u on the polar cap and associate a moving ion zone with 
the formation of a subpulse. Then it is possible to envisage 
organized circular motion in which an ion-emission zone of 
angular width 54'i{u) is followed by a proton-emission zone 
of K54>i(u). Thus at any instant, 

^ 54>r{u) (l + k,{u)) = 27r, (23) 

i 

with i = l....n in which n is here the number of elements on 
the carousel and each K is determined according to equation 
(19) by the preceding ion-zone. Comparison with measured 
values of the longitudinal subpulse separation P2 indicates 
that the order of magnitude of model values should not be 
large, Ki ~ 3, but this order of magnitude is associated 
with values of the acceleration potential, on the ion-zone flux 
lines, that are unlikely to enable significant electron-positron 
pair creation. This is then consistent with mechanism (b), 
based on protons and ions, for the production of coherent 
radio emission. 

It is easy to see that a steady uniform state of equation 
(23) with 54> = Tp4> and constant 4> that is independent of 
u does not, in general, exist. However, we have to bear in 
mind the limitations of our model based on equation (16), 
which is local in u, and on its very elementary approxima- 
tion given by equation (17). Thus in a model with a more 
realistic diffusion function fp and without the finite element 
structure, we would expect that at constant (j), K would in- 
crease as a function of u, starting from a negligible value at 
Wo and reaching a maximum before declining as a proton- 
emission zone is entered. The presence of this extremum is 
a possible basis for the existence of short-term quasi-steady- 
state solutions of equation (23) within a finite interval of u 
at Uc which also coincides with conditions suitable for the 
growth of the instability giving observable emission. It may 
be that the tendency to form a cluster, present in our model, 
indicates that a more realistic model would counteract the 
dispersion inherent in equation (23) and maintain subpulse 
shape, but it is not possible to assort that it would be so. 

However, it is possible to envisage this motion on a polar 
cap of any plausible geometrical shape, elliptical or approxi- 
mately semi-circular. The motion is, of course, independent 
of E X B drift and could be in either direction, as is seen 
in the survey of Weltevrede, Edwards & Stappers (2006b), 
or may even be bi-directional at some instant on different 
areas of the polar cap. There are many ways in which po- 
tential fluctuations could disrupt this organized motion. We 
have seen in the previous Section that a whole-surface tem- 
perature Ts large enough to give low values of V favours 
the occurrence of nulls through upward potential fluctua- 
tions. The most simple fluctuation is one that increases the 
potential over much of the polar cap and in particular at 
u = Uc- The functions es{V), and hence K, increase quite 
rapidly with V . There are two consequences. Firstly, the ex- 
ponent in the expression for the mode growth rate (Paper 
IV, equation 17) is dependent on 7^^^^ so that the condi- 
tions necessary for observable coherent radio emission are 
not reached and a null occurs. Secondly, the (unobserved) 
motion of the ion zones is substantially changed. 
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It is possible to see the nature of these changes by con- 
sidering how two variables, the proton density q{(f), t) at the 
top of the atmosphere and the ion current density satisfy- 
ing equations (17) - (19) on a section of the path at m = Uc, 
evolve over a sequence of times. In doing this, we assume 
that the ion-zone is able to move continuously and abandon 
the fixed finite-element calculation of the model described 
in Section 3.1 The ion-zone motion during a null shows how 
a form of null memory occurs. 

It is convenient, for this problem, to represent the den- 
sity q in units of pGj{0)cTp and in units of pgj{0)c. The 
ion-zones at time t have ~ 1, width 50^^' and move with 
velocity <^'^' = 5(j)^^^/Tp. Each generates a proton density 
q{(j>, t) at the top of the atmosphere which is depleted by 
the proton-zone current density = 1. Thus at an instant 
t, and as a function of 0, the proton density following an 
ion zone rises linearly to a maximum of -^f^' - 1 in a dis- 
tance Scj)'^^^ and falls linearly to zero in a further distance 
- 1)J^(^'. The velocity of an ion zone is therefore 
determined by the gradient dq/dtj) of the preceding pro- 
ton zone. This is the quasi-steady-state motion at time t, 
at which instant, an upward potential fluctuation increases 
proton production to A''-^' . The increase in ion Lorentz fac- 
tor reduces the quasi-longitudinal mode growth rate to a 
sub-critical value and so initiates the null. A new maximum 
proton density q = K^'^^ — 1 is reached at a later time t + 2Tp 
but the system velocity remains <j>^^'> until t + K^^^Tp. The 
width of an ion zone is then reduced. In detail, and with 
reference to equations (17) - (19), this happens because the 
velocity of its leading edge is reduced whilst that of the rear 
edge remains at 0^^' until t + (K^^^ + VjTp at which time the 

ion-zone width and velocity are both reduced by a factor 

of (;?(2) _ ^(1) ly ^j^g 

ion-zone velocity then tends to 
its new steady-state value = 0<^'(1 + ^'^^)/(l + K'-'^^). 
Hence subpulse drift does not stop in the model, but is much 
reduced after a null time interval of K^^^Tp. Thus for short 
nulls, there is a subpulse memory, but it does not corre- 
spond precisely with the conservation of subpulse longitude 
over a short null that has been observed, for example in PSR 
B0809-I-74, by van Leeuwen et al (2002). 

We arc not unduly disturbed by this, because there arc 
many possible spatial distributions of potential fluctuation 
over the polar cap other than the uniform case considered 
above. It is also not clear that the precise form of mem- 
ory seen in B0809-I-74 is universally observed. The upward 
potential fluctuation that caused this could reverse at any 
time during or after the sequence described in the previous 
paragraph. In relation to null memory, the question then is, 
when does observable coherent radio emission recommence? 
The problem arises owing to the very elementary nature of 
our model dcflncd by equations (17) - (19) and is that Scf)^'^'' 
in the model can be very small. The acceleration potential 
within a very small interval of (f> may be too large to permit 
adequate growth of the mode and the conditions discussed 
briefly in Section 4.1 may not be satisfied. 

5.2 Subpulse drift 

Published observations on mode-changes and subpulse drift 
demonstrate the extremely heterogeneous nature of these 
phenomena. Some features are quite distinct in a small num- 
ber of pulsars, but much less distinct or not present at all 



in others. For this reason, we have used the results of the 
extensive survey made by Wcltevrede ct al (2006b) based on 
observations of 187 pulsars selected only by signal-to-noise 
ratio, and list some of their conclusions below. 

(i) Roughly equal numbers of pulsars have subpulse drift 
to smaller or greater longitudes, and direction reversal is 
observed. 

(ii) Bi-directional drifting is observed in a small num- 
ber of pulsars. J0815-t-09 and B1839-04 have mirrored drift 
bands. 

(iii) The drift rate can be mode-dependent, having one 
of a number of discrete values of the band separation P3, 
and nulls may also be confined to a particular drift-mode. 

(iv) Drift bands can be curved or non-linear with pulse 
longitude-dependent spacing of subpulses: they are often in- 
distinct and can be found only by two-dimensional spectral 
analysis. 

(v) The band separation P3 is uncorrelated with age 
(P/2P), or with P and B. 

Conclusions (i) and (ii) are quite naturally consistent 
with the model. We emphasize again that E x B drift is 
not involved and there is no reason why, given the chaotic 
state of the polar cap, there should not be reversals or even 
bi-directionality. The nature of the chaotic state also means 
that it is unrealistic to suppose that these phenomena can 
be predictable. Conclusion (iii) is also a natural consequence 
of medium time-scale instability in the value of the surface 
atomic number Zg and hence in the potential as mentioned 
in Section 4.2. The existence of a polar-cap area emitting 
nuclei with Zs too small to produce a significant reverse- 
electron flux is obviously associated with an upward dis- 
placement of V for times very broadly of the order of Tri- 
This is quite consistent with the observed time intervals of 
10^~* s for mode-changes. Non-linear drift bands (iv) are 
simply a consequence of non-uniformity of Ki values in equa- 
tion (23). Band separation Pj was discussed in Paper II. It is 
given by P3 = {K + l)rp. In this expression, Tp is dependent 
only on the properties of the LTE atmosphere at the polar- 
cap surface, but K may have some dependence on rotation 
period and magnetic flux density; also on whole-surface tem- 
perature Ts and hence pulsar age. It was noted in Paper II 
that the distribution of P3 is quite compact, most of the 
values listed by Weltevrede et al (2006b) being within the 
interval 1 < P3 < 10 s. The fact that the carousel path is 
restricted to being near the polar-cap periphery uo may well 
act as a constraint on the values of K and so P3 that occur. 

Finally, Kloumann & Rankin (2010) have observed ran- 
domly distributed pseudo-nulls in B1944+17 of length less 
than 7P with weak emission, but well above noise. They in- 
terpreted them as a state of no subpulse on that band of 
polar-cap flux lines from which photons can enter the line of 
sight. These occur naturally in the model. The final column 
of Table 1 gives pvn, the probability that the whole polar cap 
is in the proton phase at any instant. There would then be 
no possible way of producing particle beams satisfying con- 
ditions (a) or (b) of Section 2.3 and capable of generating 
coherent radio emission. This contributes to the total null 
fraction that is observed. The probability Psn > Pvn that 
there will be no ion-zones merely in the part of the polar 
cap which is observable along the line of sight is therefore a 
natural parameter of the model. 
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6 CONCLUSIONS 

Papers I - IV, followed by the present paper, have attempted 
to predict some of the consequences of there being two popu- 
lations of isolated pulsars having opposite spin directions. In 
many important respects, the physical state of the polar-cap 
has been found to be dependent on the sign of the Goldreich- 
Julian charge density. In the fl ■ B < case considered in 
these papers, SCLF boundary conditions have been assumed 
rather than the E • B 7^ condition of the original Ruder- 
man & Sutherland (1975) model. Then the state of the polar 
cap is unstable on time-scales of the order of Tp and of the 
longer polar-cap ablation time th- These are both many or- 
ders of magnitude longer than the characteristic polar-cap 
timc-scalc of uo/c and there can be little doubt that the 
system should be able to allow the movements of charge 
necessary to maintain the SCLF boundary conditions for 
them. We therefore need to look at published observational 
data to see if there is any real evidence for the existence of 
two populations. 

The Weltevrede et al (2006b) survey is an obvious start- 
ing point. It lists 187 pulsars selected only by signal-to-noise 
ratio and examines both subpulsc modulation (the wide vari- 
ations of intensity at a fixed longitude in a sequence of ob- 
served pulses) and subpulse drift. Subpulse modulation is an 
almost universal characteristic and we assume it to be a di- 
rect consequence of the plasma turbulence that itself is now 
widely believed to be the source of the radio cmisssion (sec 
Asseo & Porzio 2006). Subpulse drift with measurable val- 
ues of both longitudinal separation of successive subpulses 
P2 and band separation P3 was detected in 72 pulsars and 
these can be compared with a set of 113 which do not show 
detectable subpulsc drift. The distributions of both sets as a 
function of age are wide and are broadly similar except that 
only 7 of the 72 have an age less than 1 Myr as opposed to 
31 of the set of 113. The distributions as functions of the 
parameter X whose critical value = 6.5 defines the CR 
pair creation threshold are also wide but the set of 113 has 
an excess &l X > Xc- There are 35 pulsars with X > Xc as 
opposed to only 7 of the 72. This is consistent with, but does 
not prove, the existence of two populations which separate 
as X moves with age to values below Xc- We propose that 
the n • B < set show the phenomena of mode-changes, 
nulls and subpulse drift as they age through growth of the 
ion-proton beam quasi-longitudinal mode (mechanism (b) of 
Section 2.3 but sec also the comments there on more obscure 
processes of pair creation for pulsars of this spin direction). 
Pair creation by ICS photons in the f2 • B > case has been 
very fully investigated by Hibschman & Arons (2001) and 
by Harding & Muslimov (2002, 2011) and our assumption is 
that this is the source of coherent emission in the population 
that has no nulls or subpulse drift. 

Having divided the Weltevrede et al list into two popu- 
lations, we should compare these with the extensive survey 
of nulls made by Wang et al (2007). No pulsar in Table 1 
of Wang et al appears in the earlier Weltevrede et al list, 
presumably because they were not known or did not sat- 
isfy the selection criteria. Of the 46 pulsars with measured 
null fractions given in Table 2 of Wang et al, 18 also do not 
appear in Weltevrede et al, but 20 are listed with a mea- 
sured value of P3 and 8 have no detected P3. A smaller but 
more recent list of nulling pulsars in the paper of Gajjar, 



Joshi & Kramer (2012) includes a further 8 that were not 
considered by Wang et al, of which 2 have measured values 
of Pa but 6 do not appear in the paper of of Weltevrede et 
al. Of the 8 pulsars that have nulls but no detected P3, 4 
have null fractions 1 — fon < 0.01 and B0656+14 is the spe- 
cial case described in Section 4.1. There are only 3 pulsars 
having substantial null fractions but no detected P3. These 
are the otherwise unremarkable B1112-I-50, B2315-I-21 and 
B2327-20. But it is not necessarily correct to regard them as 
anomalous because, as noted by Weltevrede et al, drift bands 
are often indistinct and detectable only by two-dimensional 
spectral analysis. We suggest that the results of this com- 
parison arc by no means inconsistent with our division of 
the Weltevrede et al pulsars into two populations. 

A different, but perhaps more anomalous set are the 
three pulsars known to exhibit nulls of very long duration, of 
the order of 10 d. The first of these to be found (B193H-24; 
Kramer et al 2006) enabled spin-down rate measurements 
to be made separately in both on and off states of emission. 
This was followed by similar measurements for J1832+0029 
(Lorimer et al 2012) and J1841-0500 (Camilo et al 2012). In 
each case, the off-state spin-down rate was about half that 
of the on state. But these pulsars have quite large values 
X = 3.6, 2.5 and 6.6 respectively, close to the critical value 
Xc = 6.5, and may quite possibly support self-sustaining CR 
pair creation during the on but not the off state. The conse- 
quent difference in the flux and nature of particles passing 
through the light cylinder appears to be the only physically 
plausible mechanism for a spin-down torque change of this 
magnitude. But the 10® s time-scale for both on and off 
states of emission appears a little too long compared with 
times of the order of lOr^i found from equation (22). The 
question of quasi-periodicity with such time-scales should 
also be addressed, as in the case of the RRATs (Palliyaguru 
et al 2011). We have stressed that our model is not ran- 
dom, but deterministic. Quasi-periodicitcs arc therefore not 
impossible in principle, but remain quite difficult to explain. 

Changes in the polar-cap acceleration potential have 
little direct effect on the spin-down torque. Changes in the 
proton-ion composition occurring in the process (b) nulls 
described in section 4.2 would lead to a change in the mean 
charge to mass ratio of particles crossing the light cylin- 
der, but it must be doubted whether mechanism (b) could 
explain the large observed difference in spin-down torque. 
However, it is not yet clear that this is a universal feature of 
nulls. The only other pulsar for which measurements have 
been made (PSR B0823-F26: Young et al 2012) has a frac- 
tional upper limit of 0.06 for the change in spin-down torque. 
It has a value X = 2.6, rather below the putative critical 
value Xc ~ 6.5 denoting the pair creation threshold, so that 
the relatively small change in torque could be consistent ei- 
ther with relatively weak pair production or with process 
(b). 

In the absence of electron and positron densities large 
enough to adjust and so cancel an electric field, as in mech- 
anism (a), we have to remember that acceleration or de- 
celeration remains at higher altitudes beyond ~ 10 as a 
consequence of natural flux-line curvature. But at this stage, 
growth of the quasi-longitudinal mode has already occurred 
so that further acceleration would have negligible effect on 
the coherent radio emission. This is a further distinction be- 
tween the two populations. In the f2 • B < case, the polar 
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cap of open magnetic flux lines may well have an approxi- 
mately semi-circular shape. In this case, the observed drift: 
of subpulses would be along a diameter instead of an arc of 
a circle. 

Although all the above problems remain, the model 
f2-B < polar cap described here has some positive features. 
It does not require that neutron-star magnetic fields lie in a 
particular interval. This is important because radio pulsar 
inferred polar fields can vary by up to six orders of magni- 
tude. The physical processes in electromagnetic shower de- 
velopment or in the photoelectric transitions exist in the 
zero-field limit and do not change in any qualitative way 
with increasing field. 

The model is deterministic, but chaotic, but is incom- 
plete in that its use of finite elements for the reason stated 
at the end of Section 3.1 and the very elementary nature of 
the approximation made in equation (17) appear to preclude 
the spontaneous appearance of subpulse drift from the ran- 
dom initial state which is used. It is possible to do no more 
than assert that the model can support subpulse drift in a 
quasi-stable way. 

It is unfortunate that quantitative model predictions 
depend so much on surface atomic rmmbcr and particu- 
larly on whole-surface temperature Ts, parameters which 
are not well-known. Cooling calculations (see the review of 
Yakovlev & Pethick 2004) show that Ts falls steeply at ages 
greater than 1 Myr. Bearing in mind that, for a neutron- 
star with the mass and radius assumed here, the observer- 
frame temperature is « O.STg, we can see that the 
temperatures assumed in Table 1 fall well below the val- 
ues currently observable. Although cooling in this interval 
is photon-dominated, the whole-surface temperature must 
be regarded as very uncertain. 

But having made these reservations, the model does 
represent a physically-realistic framework for understanding 
the reasons why RRATs and the varied phenomena of mode- 
changes, nulls and subpulse drift appear during neutron-star 
aging. It may be that further observations at frequencies 
below 100 MHz will provide evidence for the existence of a 
population emitting by process (b) and having spectra bi- 
assed towards lower frequencies than those for the process 
(a) population. 
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